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O I- Abstract. A new revision of the gamma flux that we expect to detect in Imaging Atmospheric 

D 1 Cherenkov Telescopes (IACTs) from SUSY dark matter annihilation in the Draco dSph is presented 

using the dark matter density profiles compatible with the latest observations. This revision takes 
ly-j . also into account the important effect of the Point Spread Function (PSF) of the Cherenkov tele- 

■ scope. We show that this effect is crucial in the way we will observe and interpret a possible signal 

profile in the telescope. Given these new considerations, some light can be shed on the recent signal 
excess reported by the CACTUS experiment. 
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Nowadays, the indirect search for supersymmetric (SUSY) dark matter is possible by 
means of the new Imaging Atmospheric Cherenkov Telescopes (IACTs). This search is 
based on the detectability of gamma rays coming from the annihilation of SUSY dark 
matter particles that takes place in those places in the Universe where the dark matter 
density is high enough. IACTs in operation like MAGIC or HESS, or in the near future 
the GLAST satellite or GAW, will play a very important role in this kind of SUSY dark 
matter searches. 

The expected total number of continuum y-ray photons received per unit time and 
per unit area, from a circular aperture on the sky of width a t (which represents the 
resolution of the telescope) observing at a given direction relative to the centre of the 
dark matter halo is given by: 

F(E>E th ) = ^-f SUS Y-U^o), (1) 



fsusY = u(?o) = J J{y)B(a)dn, 
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where the factor fsuSY encloses all the particle physics, and the factor U (*Po) involves 
all the astrophysical properties (such as the dark matter distribution and geometry con- 
siderations). This astrophysical factor also accounts for the beam smearing, where: 



f pj m (r)dl, dl = ±rdr/Jr 2 -dlsm 2 y (2) 

Jl.o.s * 

is the integral of the line-of-sight of the square of the dark matter density along the 
direction x ¥, and B(Q.)dCi is the Gaussian beam of the telescope: 



B(Q.)dQ. = exp 



sinQdOdq. (3) 
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This last factor, commonly known in the astrophysical community as the Point Spread 
Function (PSF) of the instrument, plays a very important role in the way we will "see" 
a possible signal in the telescope. However, most of previous work in the literature did 
not take into account its effect (except O]). 

A very important question concerning the indirect search of SUSY dark matter is 
where to search. Because of the fact that the factor U (^o) gamma flux is proportional 
to the square of the dark matter density, we will need to point our IACTs telescopes 
to places with a high concentration of dark matter. In principle, the best option seems 
to be the Galactic Center (GC), since it satisfies this condition and also it is very near 
compared to other potential targets. However, the GC is a very crowded region, which 
makes it difficult to discriminate between a possible y-ray signal due to dark matter 
annihilation and other astrophysical sources. 

There are also other possible targets with high dark matter density in relative proxim- 
ity from us, which are not plagued by the problem of the GC, e.g. the Andromeda galaxy, 
the dwarf spheroidal (dSph) galaxies - most of them satellites of the Milky Way- or even 
huge cluster of galaxies (e.g. Virgo). dSph galaxies represent a good option, since they 
show very high mass to light ratios, and at least six of them are nearer than 100 kpc 
from the GC (Draco, LMC, SMC, CMa, UMi and Sagittarius). Draco is the dSph with 
strongest observational constraints, located at 80 kpc. This fact is very important if we 
really want to make a realistic prediction of the expected y-ray flux. Moreover, recently 
the CACTUS collaboration reported a possible gamma-ray excess from Draco, which 
makes a detailed study of this dSph even more attractive and necessary. 



DARK MATTER DISTRIBUTION IN DRACO 

In our modelling of Draco we used the sample of 207 Draco stars with measured line-of- 
sight velocities originally considered as members by 03]. We used a rigorous method of 
removal of possible interlopers originally proposed by [1] and applied to galaxy clusters. 
The method relies on calculating the maximum velocity available to the members of 
the galaxy assuming that they are on circular orbits or infalling into the structure. The 
method was shown to be the most efficient among many others for interloper removal 
recently tested on cluster-size simulated dark matter haloes by 10]. Its applicability and 
efficiency in the case of dSph galaxies was demonstrated by flj]. 



T 




1 - \ 

_J I 1 L_ 

0.1 1 10 

r [kpc] 

FIGURE 1. The best-fitting dark matter density profiles for Draco with a cusp (solid line) and a core 
(dashed line). 

TABLE 1. Best-fitting parameters for the dark matter 
profiles with a cusp (a = 1) and a core (a — 0). The last 
column gives the goodness of fit measure % 2 /N. 
profile M D /Ms rb/fl s P ]f/N 

cusp 830 7.0 -0.1 8.8/9 

core 185 1.4 0.06 9.5/9 



We assumed that the dark matter distribution in Draco can be approximated by the 
formula: 

p d (r)=Cr-«exp(--0 (4) 

proposed by [2], which was found to fit the density distribution of a simulated dwarf 
dark matter halo stripped during its evolution in the potential of a giant galaxy. [2] found 
that the halo, which initially had a NFW distribution, preserves the cusp in the inner part 
(so that a = 1 fits the final remnant very well) but develops an exponential cut-off in the 
outer parts. Here we will consider two cases, the profile with a cusp a = 1 and a core 
a = 0. It remains to be investigated what scenarios could lead to such core profiles. 

We modelled the velocity dispersion and kurtosis profiles calculated from the data 
solving the Jeans equations with dark matter profiles given by Eq.(4). The best-fitting 
solutions in the case of a cusp and a core profile are given in Table 1 . M^/Mv, is the ratio 
of the total dark matter mass to total stellar mass, r^/Rs the break radius of Eq.@ in 
units of the Sersic radius of the stars and /3 the anisotropy parameter of the stellar orbits. 

Fig-dshows the best-fitting dark matter density profiles in the case of the cusp (solid 
line) and the core (dashed line). As we can see, both density profiles are similar up to 
about 1 kpc, where they are constrained by the data. The reason for very different values 
of the break radius in both cases is the following. The kurtosis is sensitive mainly to 
anisotropy and it forces /3 to be close to zero in both cases. However, to reproduce the 



velocity dispersion profile with /3 ~ the density profile has to be steep enough. In the 
case of the core it means that the exponential cut-off has to occur for rather low radii, 
which is what we see in the fit. The cusp profile does not need to steepen the profile so 
much so it is much more extended and its total mass is much larger. 

THE ROLE OF THE PSF 

In order to compute the expected gamma flux, we need to calculate the value of the 
astrophysical factor £/(y/b), given in Eq.©, for the core (a = 0) and cuspy (a = 1) 
density profiles as given by Eq.© with the parameters listed in Table [T] The results are 
plotted in Figure El where we used a PSF of 0. 1° (which is the typical PSF for an IACT 
like MAGIC or HESS). As we can see, it is possible to distinguish between both density 
profiles thanks to a different and characteristic shape in each case. 
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FIGURE 2. Draco flux predictions for the core (dashed line) and cusp (solid line) density profiles. A 
PSF of 0.1° was used 

To illustrate the PSF effect on the shape of the flux profile measured with our IACT, 
in Figure |3] we show the same as in Figure |3 but here for a PSF= 1° (the PSF of the 
CACTUS experiment). It is clear that, although we use different dark matter density 
profiles, a worse PSF makes both resultant flux profiles indistinguishable. According to 
these results, however, one may think that we could distinguish them by means of the 
value of the flux near the centre, i.e. for small values of y/Q. However, this is not possible 
in practice due to the large uncertainties, which come mainly from the particle physics 
(uncertainties around three or four orders of magnitude in most cases, see 0]). 

For the same dark matter density profile, a worse PSF flattens the flux profile. It can be 
clearly seen in Figured where we plot the flux predictions for the cuspy density profile, 
and for two different values of the PSF (0. 1° and 1°). It is also possible to see the same 
effect if we make a comparison between Figures El and |3] 
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FIGURE 3. Same as Figure|2]but for a PSF= 1°. 
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FIGURE 4. Draco flux predictions for the cusp density profile, and for two different PSFs. Solid line 
corresponds to PSF=0.1° and dashed line to PSF=1°. 



A DISCUSSION OF CACTUS RESULTS 

CACTUS is a ground based y-ray telescope located in California. The experiment is 
sensitive to y-rays above 50 GeV and it has an effective area of about 50000 m 2 . It 
was first designed for solar observations and not for y-ray astronomy. Because of this 
fact, it has a poor PSF of around 1°. Recently, the CACTUS collaboration reported a 
y-ray excess from Draco ([5]). In [8], the CACTUS data were superimposed on different 
flux profiles, each of them related to different models of the dark matter density profile. 
However, we must note here that their flux estimations were computed without taking 
into account the important role of the PSF. 

As we can see in the left panel of Figure |51 if we adequately introduce this effect, 
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FIGURE 5. Flux profiles for the four models of dark matter density profile used in 1 8], computed using 
the CACTUS PSF= 1° (left panel), and using an improved PSF= 0. 1° (right panel). 



it will be impossible to discriminate between different flux profiles, i.e. between the 
four models for the density profile described in BSD, using the CACTUS PSF. Only the 
absolute flux could give us a clue, but as it was mentioned before, there are too many 
uncertainties in the y axis. In the right panel of Figure|5]the same exercise was done, but 
with an improved PSF= 0.1° (e.g. the MAGIC PSF). In this case we could distinguish 
between different flux profiles, i.e. different models for the dark matter density profile. 

DRACO AND GAW 

GAW is a R&D path-finder experiment, under development, for wide field y-ray astron- 
omy. GAW will operate above 0.7 TeV and will have a PSF ~ 0.2°. It will consist of 
three identical telescopes working in stereoscopic mode (80m side). The mean goal of 
GAW is to test the feasibility of a new generation of IACTs, which join high sensitivity 
with large field of view. GAW is planned to be located at Calar Alto Observatory (Spain) 
and a first part of the array should be completed and operative within winter 2008. GAW 
is a collaborative effort of Research Institutes in Italy, Portugal, and Spain. 

It is possible to make some calculations concerning the possibility to observe a 
y-ray excess in the direction of Draco by GAW. If we are only interested in flux 
detectability (i.e. no discrimination between different dark matter density profiles), we 
should calculate simply the integrated flux that we will observe in GAW due to dark 
matter annihilation in Draco using Eq.® with E t h = 0.7 TeV. To make the calculations, 
we suppose Draco to be around 1 .5° in the sky, a GAW PSF= 0.2° and aS/N > 5. There 
are, however, large uncertainties (e.g. the exact value of fsusY for 700 GeV, the Draco 
inner profile, etc). If we take fsusY = 10~ 34 ph ■ GeV 2 ■ cm ■ s~ l for 0.7 TeV (as in 
0]), and a density profile given by Eq.© with a = 1, we obtain for the total flux: 

Foraco = 2.4876 • 10~ u ph cm- 2 s~ l (5) 

Given this value, it may be possible to detect a y-ray excess from Draco by GAW, 
since the minimum detectable flux in 50 hours at 5<7 level with this IACT is 3.5 • 



1(T 12 ph cm- 2 s~ l (LJ). 



CONCLUSIONS 

We can summarize the main conclusions as follows: 

• The PSF of the instrument is crucial to estimate correctly the expected gamma flux 
profile due to dark matter annihilation in IACTs. 

• The effect of the PSF could make it impossible to discriminate between different 
models of the dark matter density profile. 

• Concerning the y-ray excess reported by CACTUS in the direction of Draco, their 
results (if real) should be interpreted carefully. There is no possibility to say, given 
the poor PSF of the experiment, if the dark matter density profile is cusp or core. 
May be they only detected an excess signal, which has to be confirmed. 

• Moreover, if the CACTUS excess is real, MAGIC should see the signal without 
problems and could distinguish among the different dark matter density profiles. 

There is no doubt that GLAST, with a PSF< 0.1°, will be very important in the 
indirect search of dark matter. In that sense, also IACTs with large field of view and 
high sensitivity are the next step in this search. The R&D experiment GAW represents a 
first attempt in that direction. 
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